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TECHNICAL NOTE 3667 


WING-BODY COMBINATIONS WITH CERTAIN GE (METRIC 
RESTRAINTS HAVING LOW ZERO-LIFT WAVE DRAG 
AT LOW SUPERSONIC MACH NUMBERS 
By Harvard Lomax 


SUMMARY 


Several variational problems involving optimum wing and body combina- 
tions having minimum "wave drags for different kinds of geometrical restraints 
are analyzed. Particular attention is paid to the effect on the wave drag 
of shortening the fuselage and, for slender axially symmetric bodies, the 
effect of fixing the fuselage diameter at- several points or even of fixing 
whole portions of its shape . 


INTRODUCTION 


Recently several authors have used linearized theory to study the 
wave drag of wing -body combinations traveling at supersonic speeds (see, 
e.g., refs. 1 to 5 ). These studies have clearly demonstrated the impor- 
tance of finding the rove drag of a whole airplane rather than the separate 
wave drags of its various parts (wings , fuselages, etc.), since the magni- 
tude of the interference terms can predominate. In effect, this means 
that various optimization problems for bodies - such as the problem of 
finding the body shape having a minimum wave drag for a given volume - 
should be re-e xamin ed when interfering wings or other bodies are in the 
same flow field. In many cases the solution to the new problem differs 
from the body-alone problem only in interpretation. 

The purpose of this report is to study minimum wave- drag combina- 
tions which satisfy a few of the many possible geometric restraints per- 
tinent to the interests of airplane designers . An attempt has been made 
to anal yze the various problems in a unified manner so that extensions 
to other kinds of restraints can be deduced. 
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LIST OF IMPORTANT SYMBOLS 
aspect ratio 

source distribution equivalent to wing in sense defined by 
equation (3) 

multipole distribution of order n 
wave drag 

portion of drag due to all tbe nth order multipoles for n>0 
See equation (8). 

additional drag resulting from restraint (See eq. (ll).) 
restraints defined in equations (19) 

distance between apexes on x axis of forecone and aftercone 
enclosing wing (See sketch (c).) 

length of basic body 

length of modification to basic body 

Mach number 

TT 2 

p o^o 

2 

average body radius 

fuselage area in cross section normal to the free stream 

normal (to free stream) projection of wing area in section cut 
by plane X]_ = x + p y 1 cos 6 (See sketch (b).) 

speed of free stream 
volume 

Cartesian coordinate system (See sketch (a).) 

source distribution representing the fuselage modifications 
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n/m 2 - 1 

polar coordinate (See sketch (a).) 

free -stream density 
See equation ( 17 )- 


velocity potential 


BASIC THEORY AKD ASSUMPTION 


Basic Theory 


Many of the discussions and derivations contained in the following 
are carried out on the assumption that the reader is familiar with the 
concepts presented in reference k which should he considered as a first 
part to this report- In particular, 1 Z 

an acquaintance with the solutions 
to the wave equation referred to as 
"multipoles" is assumed, together 
with Hayes ’ invariance principle and 
the consequent multipole distribu- 
tions equivalent to a wing in the 
sense that both induce the same 
momentum spectrum at infi n ity. 



The entire analysis used herein 
is based on the assumptions and ideali- 
zations necessary to develop the 
linearized equation for the velocity 
potential, 9 , in supersonic flow, 
namely 

P 2( Pxx “ <Pyy “ <Pzz = 0 C 1 ) 

where 0 2 = M 2 - 1 and the reference 
coordinate system 1 is shown in 
sketch (a) . The analysis is further 
restricted to the solution of prob- 
lems involving a given uncambered 


1 It should be stressed that the 

stream direction (wind axes) . 



Sketch (a 


axis is parallel to the free- 
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and untwisted wing mounted „ centrally on a vertically symmetrical fuse- 
lage, the entire configuration being at zero angle of attack. 


Additional Assumptions 


We now make the two additional assumptions: one, the value of (3 A, 
where A is the w ing aspect ratio, is small; and two, the value of 0R/L e , 
where R is the average body radius and 2Le is the distance along the 
x axis in which the multipole strengths differ significantly from zero, 
is small. 

One can evaluate the significance of these assumptions by studying 
their implications relative to the source and multipole distributions 
used to simulate the wing and body. Suppose, for example, a group of 
multipoles are placed along the body center line, their strengths, C n (x), 
being fixed by the condition that a circular cylinder in the vicinity of 
the body is a stream surface when the velocity field induced by these 
multipoles is combined with the velocity field induced by the source 
sheets representing the given wing. With the assumptions of small 0A 
and PR/L e mentioned above, the C n ’s, for n greater than 0, can be 
shown (see, e.g., ref. A) to have a negligible effect on the wave drag. 
Hence, all the multipoles (for n > 0) that combine with the wing to make 
a circular cylinder a stream surface and any additional multipoles (for 
n > 0) added to make the body have mild distortions from such a surface 
are negligible in evaluating the wave drag. Therefore, under the assump- 
tions mentioned above, out of all the singularities distributed along the 
body axis, it is necessary, in studying the wave drag, to consider only 
the sources. 

With the restrictions to small values of 0R/L e and mild body dis- 
tortion (see Ward, ref. 6, for a discussion of orders of magnitude), 
slender -body theory can be used to calculate the body shape , and on the 
basis of this theory one can show (see ref. k, Appendix B) that Sf(x), 
the body cross-sectional area measured normal to the free stream, is 
completely determined by the axial source distribution alone. Hence, if 
only the exposed panels of the wing are Used to calculate the C n ’s, C Q 
is negligible and the entire axial source distribution a^x) is related 
to the geometrical properties of the body by the relation 


Oxd( x ) 


Uo 


dSf 

dx 


UoSp '(x) 


( 2 ) 
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The Wing Equivalent Source Distribution and the 
Optimum Cancellation Sources 

Let the given wing lie in the z± = 0 plane . According to Hayes * 
theorem (ref. 7), the wing equivalent source distribution [ao(x)] w is 
obtained by accumulating on the Xi axis, at a distance x from the 
origin, al 3 the wing sources intercepted by the line xi = x + ffy^eos 0, 
and then, for a fixed x, averaging these values as 0 varies from 0 
to 2it. Thus, using thin-airfoil theory to relate the planar source sheet 
to wing geometry, one finds 




( 3 ) 


where Sv/(x,0) = d/dxts^x,©) ] and 
Sw(x,0) is the normal (to the x axis) 
projection of the wing cross-sectional 
area intercepted by the plane 2 
x x = x + fryicos 0 as shorn in 
sketch (b) . Without the addition of 
further restraints, the optimum source 
distribution along the x x axis is 
that which just cancels the wing- 
equivalent source distribution. Fur- 
ther, this can be interpreted directly 
in terms of both fuselage end wing 
geometry by means of equations (2) yj 

and (3)* Thus, with no further 
restraints, the best fuselage shaping, 
for a wing-body combination satisfying 
the assumptions discussed above, satis- 
fies the equation 

i r 25t 

Sf ’(x) = - ^ J Sw’(x,0)d0 00 

and h&s any reasonably smooth cross- 
sectional form. Notice that the total 


S(x t $ ) - Normal projection 
of wing area along AA 




Sketch (b) 


2 The true oblique plane is given by the equation 
x x = x + pyjcos 0 + pz x sin 0 

but, to be consistent with the assumptions basic to linearized theory, 
the variation with z x is neglected. 
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volume taken out of the fuselage is exactly equal to the total volume of 
the exposed portion of the wing. Hence , the total volume of the modified 
comb in ation is the same as that of the original smooth cylinder. 


The Drag 


The total wave drag of a system can be expressed in terms of its 
actual or equivalent multipole distributions as 

00 

D = 2D 0 + ^ Dn (5) 

1 

where Dq is the drag caused by the nth order multipoles a n (x) and is 
given by the equation 


Ih 

q. 


jSD. 




JJQ - LJ 

L 


<3x2an^ n+1 ^ (xi)a n ^ n+1 ^ (x 2 ) In | Xi - x 2 | , n = 0,1,2 , 


( 6 ) 


where a n ^ n+1 ^(x) represents (d/dx) a n (x) . Under the. assumptions 


given above, the magnitude of E Dq is small . Let us designate it by 
D e , so that, in general. 


D = 2Dn + Dc 


(7) 


On the other hand, the total wave drag of a system composed of the 
combination of a wing and a body can also be written symbolically as 


D = Dw + ^wb + 


( 8 ) 


where 

Dw = drag of the wing alone 
= drag of the body alone 
= interference drag 



MCA TN 3667 


7 


The various components of wave drag defined, in equations ( 7 ) and (8) 
help one to evaluate more readily the drag reductions that can be realized 
from appropriate fuselage formations. Thus, if the fuselage shape satis- 
fies equation (^) , the total wave drag of the combination under the assump- 
tions that f3A and 0 R/l 6 are small can be written either as 


D = D e 


t9) 


or as 


D — - Dfo 


( 10 ) 


If, in finding the fuselage shape, 

(a) the multipoles representing a wing and a body flying separately 
are assumed to represent the same wing and body when combined (i.e., 
the shape fields can be superimposed) , 

(b) the multipoles representing the fuselage are equal in magnitude 
but opposite in sign to the wing equivalent multipoles, 

then equation (10) holds without the assumption of small fiA and 0R/L e . 

In subsequent problems we will discuss the effects on the wave drag 
and fuselage area distribution of adding certain additional restraints to 
the body geometry. The addition of such restraints may or may not change 
the relation given by equation (10) , but they must always add a term to 
equation ( 9 ) so that 


D = Dg + Dre 
Dre > 9 


(H) 


WINGS CENTRALLY MOUNTED ON SLENDER QUASI -CYLINDERS 


This section is devoted to the solution of two problems involving a 
given uncambered and untwisted wing mounted centrally at zero angle of 
attack on a tube that is cylindrical forward of some point ahead of the 
w ing . The problems are, in both cases, to find the area distribution of 
the fuselage behind the cylindrical portion that will minimize the wave 
drag of the combination. 
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Shortening the Fuselage 


Remembering the assumptions listed at the beginning of this section, 
let us consider the following problem: 

(i) Given a wing and a slender fuselage having the same normal 
cross-sectional area in all planes ahead of the plane x = -L 0 ' 

(see sketch (c)), what is the optimum fuselage area distribution 
behind the plane x = -L 0 ' if the fuselage must end at the plane 
x = l 0 ? 



where from the conditions stated in the problem and the geometric inter- 
pretation to the fuselage sources given by equation (2), <x 0 (x) is zero 
for values of x outside the interval 3 -L n ' < x < l Q . 

3 It is necessary, for equation (6) to be valid, that (^(-Lq 1 ) and 
GCo(l 0 ) he zero. This implies that clq(x) must be continuous and if the 
body shape is given by equation (2), this, in turn, means that the stream- 
wise gradient of body cross-sectional area must be continuous. It was 
pointed out in reference k that a 0 (-L 0 l ) and a 0 (L 0 ) will both be zero if 
the wing has no blunt edges along which the normal component of the free- 
stream Mach number is unity or greater. 
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Consider now a variation of equation (12) for a fixed a Q (x) in the 
interval -L 0 ' < x < L 0 and a free variation of Oo(x) in the sub- 

interval -L 0 * < x < l 0 . There results 



2 r , r° . 

J <3xi Ba o(xi) J t a 0 1 (x 2 ) In I Xx - x 2 | dx 2 

0 -L~ ' “Ln 


L o 

/ ’ (x 2 ) ln I Xx - x 2 1 dx 2 

-L ' 


= 0 


Integrate once by parts with respect to Xx (since the variations bct^-Lo') 
and Sa 0 (Z 0 ) must be zero) . Then, by the fundamental lemma of the calculus 
of variations, the bracketed term must be zero for -L 0 ’ — Xl — ^o and one 
finds the condition 



a 0 , (x 2 )dx 2 
xi - x 2 x 



a 0 ’(x 2 )dx 2 
Xx- x 2 


-L 0 '<xi<Z 0 (13) 


Equation (13 ) is an integral equation which can be inverted (by methods 
such as those outlined in ref. 8). Inverting, integrating once with 
respect to x, and expressing Oo(x) and a Q (x) by means of equations (2) 
and (3), respectively, one finds 


Sf’(x) 


1_ 

2it 



7(^o- x)( L 0 * + x) 



Sy/Qcx^e) 

(xx - x) \J (L q ’ + Xl) (xx - l Q ) 


(Ik) 


which gives the optimum fuselage area distribution under the conditions 
and assumptions posed. 

The wave drag of the combination represented by equation (l4-) can be 
expressed either in the terms defined in equation (8) or (ll) . Let -us 
first consider the form given by equation (8) . If the expression for the 
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drag of an nth order multipole distribution is integrated once by- 
parts, there results since a rL ' n+1 ^ (-Iq 1 ) = a n n+1 (-^o) = 0 


Bn 

1 


52 n 


IfrtU r 


/° 

-L ’ -L„* 


®n 


(n+i) 


(x 2 )dx s 


Xi -x 2 


Using this expression, one can readily show that equation (12) yields 


®v/b 


"o “o 

= — ^-5 / a 0 (x i )dx 1 / 

2*U * • 


o _L, 


o 0 , (x g )dx 2 
Xi - x 2 


1 

2rtU, 


l Q L o 

TS / a 0 (x 1 )dx i / 

O '-'t ' '-t ‘ 


a 0 , (x 2 )dx 2 

Xi - X 2 


so that, by equation ( 13 ) 


f a 0 (xi)dxr r 

1 2flU ° J -L , ’ -L' 


a 0 , (xs)dx 2 Djj 


Xx -x 2 


Hence, for any combination satisfying equation (l*0 , once again 


D — - Dfc> 


On the other hand, D^, the increase in drag caused by shortening 
the fuselage can also be obtained. Integrating equation (12) by parts. 
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one has (note D^g = 2D 0 ) 


Lo 


Pre 

q. 


- ^ f ta 0 (xi) + aoWJbc, f ta °' (X2) ± 
^^o J T t j r ' Xi - x 2 


Combined with equation (l3)j this becomes 


Pre 


q 2«U 0 ^ 


j->o ■*- 

a Q (xi)dxi ^ 
-L~ 


a 0 ‘(x 2 ) + cio'Cxs) 


xi - x 2 


dx 2 


The derivative of equation (l4) with respect to x gives 

Lr 


I ' (x) + Oo'(x) = — T7?^-p.r=-Z=l f 

. jin/ (Zq - x ) (L 0 * + x) 


0 a 0 l ( x i) n/(Lo ' + Xi) (xi - 2 0 ) 


X - Xi 


dxx; 


-Lq < x < 


so 


Pre 

1 “ 2jtU, 


— 

b f a o(^)<^a 




a Q '(x 2 )dx 2 
x x -x 2 


-L, 


dx 2 


, (xi-x 2 )V (L 0 ’ + x 2 )(Z 0 - x 2 ) 


/ 


a. 0 '(x 3 )JiL 0 ' +x 3 ) (x 3 -l Q ) 


<1X3 


x 2 - x 3 


which reduces to 


Dl 


re 

q. “ 2«U ( 


XJQ j- 

5 ? / a ° (xi)dXl / 


a o ’ (x 2 ) I (L 0 > + x 2 )(x 2 - z 0 ) 
(L 0 * +Xi)(xi - I 0 ) 


dx 2 ( 15 a) 


*0 x,-x 2 
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or, alternately. 


fire 

1 


-1 

2*U 0 2 



a 0 '(xi)a 0 '(x 2 ) 


In 


xi - x 2 1 (l 0 + L 0 1 ) 


[«7 (xi + L c 1 ) (x 2 - l Q ) +-J (x 2 + L 0 1 ) (xi - Z Q ) ] 


dxidx 2 (15b) 



Constrained Fuselage Areas 


Another class of problems is that 
in which the magnitude of the fuselage 
area is fixed at various points . Sup- 
pose, for example, that a fuselage 
shaped according to equation (4) had 
in some region a cross-sectional area 
too small to be acceptable for some 
practical purpose. The question is, 
then, what\is '-the best shape for given 
values of minimum fuselage cross- 
section area at given planes and what 
is the penalty in wave drag caused by 
the added constraints?- Before consid- 
ering the general case of an arbitrary 
number of restraints, let us first con- 
sider the simple problem: 


(ii) Given a wing, what (under the various assumptions given 
above) is the area distribution of the adjoining fuselage which 
has a prescribed area at three given stations (the initial, the 
final, and an intermediate station x = dij see sketch (d)) and 
yields a mini mum wave drag for the combination? 


As before, let a Q (x) represent the wing equivalent source distribu- 
tion. Then the drag caused by the restraints can be written 


fire 

<1 


1 

2*U 0 2 



[a 0 (xi) +U Q Sf 


r L ° 

(x x )]dxi / 

-L„ 


a 0 '(x 2 ) + U 0 Sf" (x 2 ) 
dx 2 


(16) 


Xx - x 2 
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where Sp(x) is the unknown fuselage area to he optimized. For simplic- 
ity, replace the unknown Sp(x) by cf(x) where 


x 2 Jt 

( x ) = J ~~~ + Sf(x) = ^(x,9)d.e + Sf(x) 

-Lo ° 0 


d(-L 0 ) - S 0 

ff (^o) = S 2 


Let 


0 o (x) = d(x) ; -L 0 < x < di 

<*i(x) = cr(x)j di < x < L q 

and the restraints on the fuselage area give the relations 

= Jn 


n M 

2 Jt 
r\ 

/ Cq ' (x) dx = 

Sx - s 0 + -Z- s^d^ejde 

-Lo 

0 


r L ° 

r 23t n 

r\ 

/ Ci'(x)dx = 

s 2 - s - — / SwCd^eJde 

y ii 

2 * J o 


> (IT) 


(18) 


(19a) 


= Ji - J 0 (191) 


where J Q and J x axe constants fixed by the given constraints. Notice 


Ji = (S 2 - S 0 ) 


(19c) 


so the constant jy is a measure of the difference between the initial 
and final areas . 


Using the usual variational techniques, we can write, 
tity to be minimized. 


d i , Lo 


for the quan- 
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or 


_£ 


/ dl 

-Lo 


cr 0 '(xi)dxi 


/ 

L-Lo 


di 


a 0 ll (x g )dx 2 
Xi -x 2 


■t 

di 


di" (x 2 )dx g 
xi - x 2 


+ \ 



(xi)dxi 


r 


d o" (x g )dx g 
Xi -x 2 



di" (x 2 )dx 2 + ^ 
xi - x 2 


By taking the variation and satisfying the conditions at the end points , 
one obtains the two s imul taneous integral equations 


r 1 dp” (x 2 )dx 2 + r 

x x -x 2 J 
-L 0 di 


di"(x g )dx 2 
xi — x 2 



-L 0 < x x < di] 


> ( 20 ) 



dp” (x 2 )dx 2 
Xi -X 2 



di" (x g )dx g 

Xi-X 2 



di < Xi < L 0 


The set of equations (20) is identical to that analyzed by Adams 
(ref. 9 ) , page 1 ^, for bodies of revolution with fixed areas at the 
initial, final, and an intermediate section. In the interest of subse 
quent generalization, however, we will consider its solution in the 
following way: First write the equations (20) in the equivalent form 



a" (x g )dx g 
xi -x 2 



-L 0 < xi < di 


di < xi < L 0 


( 21 ) 
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One can show that 


2 

« / \ A + Bx^ -j L n — cLtXo 

(xa = j L 3. -- • cicosi1 . 1.- . 1 


^ -V <= 

■o - x 2 


L 0 |x 2 - di| 


( 22 ) 


is the solution to the integral equation (where A, B, and Cx are con- 
stants) since 

r 




(x 2 )dx £ 


B - ClCOS -1 hjri 

V'O/ 


-Iv 


Xl- X 2 


-< 


-It 


B+ Cxcos 


-1 


-d 


-L 0 < xi < di 


dx < Xi < Lr 


(23) 


which satisfies equation (21) . The constraints can now be satisfied by 
means of the equations 


/ x . 

cr"(xx)dxx = A cos -1 f^-J - B J L 0 2 - 

-Ln 


P 

X + 


Cx 


(dx - x)eosh -;L -^2 ^L 

L 0 x - dx 


-s/Lq 2 - dx 2 cos" 1 ^^ 


and 


a(x) - S 0 = A 


(24) 


X cos 


+ *7l 0 2 - X3] - ]| j^x ^Lo 2 - x 2 + Lq 2 cos -1 


Cx 

2 


(dx - x) 2 cosh 


-1 Lq 2 - dl X 
L 0 |x - dx | 


\Il q 2 - dx 2 (dx - 2x)cos -1 ^^ + 


<s/(L 0 2 -di 2 ) (L 0 2 - x 2 ) 


(25) 
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Since 0*(-L o ) = 0(L O ) = 0 



Solving for Ci and B, ye find 



From equation (25) the fuselage cross-sectional area can be written 





2 It / 

3f(*) = So f BvU,e)d9 + U 


/ - 
ItJb 


(L 0 2 - d 1 x>J(L 0 2 - x 2 ) (Lo 2 - d x 2 ) - 


2 , L-,” 5- - dix 

L 0 2 (x - d x ) cosh -1 - 

L 0 x- dx J 


+ 


Ji|^(Lo £ “ dx 2 ) 2 coa - 1 (~) - *] (Lo 2 rx 2 ) (Lo 2 -di 2 ) jL 0 2 -di' 


(L 0 2 - dix)cos -1 (rp^i) + (x-di) 2 dx-v/L^ 
\- L, o / J 


-o 2 “ <*i S 


+ 


_ 2 , i Lo ” .xdi 

Lo cos ^jcosh a -|J 


} 


In terms of the wing, body, and interference drag components defined in equation (8) , 
total wave drag is 


D Dy Db . 4 

5 = f ‘ T + (s °‘ 32)b + Cl 


itSx - SqCos ~ Sacos" 


<£) ] 


where B and Cx are given above. The equation for is 


= ~77 "g"T~2V l ^o 2 Jo 2 - SitJoJil dxN/Lo 2 - dx 2 + L 0 2 cos“ 1 (lii 


^oJ J 


\ 

Ji 2 -Jlq 2 - dx 2 + 2dxVLo 2 - d'x 2 cos 1 + L o 2 c ° 3 “ 1 (x^) j ^ 


(di-x)+ 


(26) 

the 

(27a) 


(27b) 


w 
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If the additional specification is made that the initial, and f inal 
areas are the same, the solution simplifies considerably, since, for such 
cases, Ji = 0 and equations (2 6) and (27) reduce to 


/ 2Jt 

S w (x,0)d0 


(l 0 2 - d^r 


Si - So + 


2 *± ‘ 


S w (d x ,0)d0 


(L 0 2 - dxxW(L D 2 -x 2 ) (Lq 2 - d x 2 ) - L q 2 (x - d 1 ) 2 cosh-i -2 - 7 

L 0 x - d x _ 


(28) 


and 


B _ Bv/ 

' 1 q. 



2itL Q 2 

( L 0 2 -^ 2 ) 2 


(Si - S 0 ) 


Si - Sq 


+ 


1 

2« 





(29a) 


q. 


itL, 


(L 0 2 - d x ) L 


Sx - S 0 


+ 



( 29 b) 


Often the exact statement of the restraint is that S(x) shall not 
he less than Si at x = d x . In such cases care must he used in applying 
equations (26) and (27) or (28) and (29) , since they are only valid when 
the fuselage cross-sectional area at d x is exactly S x . If such Is the 
case, equations (26) and (28) give the optimum body shape only if Jo > 0, 
that is, only if 


Si > S 0 - 



Otherwise the optimum variation of area is given by equation (4) . 
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Next let us generalize the analysis leading to equations (2 6 ) and ( 27 ) 
by considering the following problem: 

(iii) Given a wing, what is the area distribution of the adjoin- 
ing fuselage which has prescribed areas at n + 1 stations (includ- 
ing the initial and final ones fixed by the Mach forecone and 
aftercone enveloping the wing, see sketch (e)) and yields a mini- 
mum wave drag for the combination? 


/\ 


By analogy with equation (22), the '•’q 

integral equation for cr" (x) (where cr(x) Modified 
is defined by eq. ( 17 )) that must be body 

satisfied for a m inim um wave drag can be /K 

witten at once in the form 


Original 

cylinder 


J 

I 


a" (x 2 )dx 2 
Xi - x 2 


= d-£ < xi < d^ +1 , 

i = 0,1, . . ., n (30) 
The 



where d Q = -L 0 and d n+1 = L 0 
quantity 


- 1 -* 


A_+ Bx 2 \'r.^^v,-l L Q ~ X 2 d i 


XI 

( 31 ) 


/ Areas S Q , 

s i ' s n , 


S n + J> and 
wing all 
’* fixed. 

Sketch (e) 


is a solution to equation (30) since it yields 


v-i 


I 4 * 

"ho 1 V 


d v-1 < Xi < d v , v = 1,2, . . ., n + 1 (32) 


n 


I '* 006-1 (D - 


A_p 

n+i 


in which 
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Further, it is apparent from equations ( 24 ) and ( 25 ) that, with the con- 
ditions 0'(-L o ) =0 , (L O ) = 0 


0’(x) 


-B «n/l o 2 - x 2 + 



-x)cosh -1 


L 0 2 - xdj 

L o|x-di| 


( 33 ) 


and 


2JT 


Sf (x) = S c 


-hi 


s„(x,e)ae - | 




+ Lr^cos -1 ! — 

Lo. 


Ci^di - x^coslT 1 L ° 2 “ Xd± 


i=i 


Lq|x -dil 


- \JL 0 2 - dj_ 2 d^cos -1 ^^+N/L 0 2 - x 2 


( 34 ) 


The wave drag due to the restraints can be obtained by using equa- 
tions (32) and (16) . Thus 


fire 

q. 


n 

= §• B (°o ~ ^n+i) + \ Yj, Ci 11(11 " a o cos ~ 1 Q~) 

i=i 


I - d n+1 cos 


-1 


-di 


J o/ J 


( 35 a) 


or in terms of the components defined in equation (8) 


D _ Dw 

q. “ q 


Db 


+ B(S 0 


>*I* 


jtSj_ 


1=1 




(35b) 


where 0 ± = 0(d ± ). Notice = S Q and 0 n+1 = S n+1 , so when S D = S^, 


^p = | - so) 

i=l 


( 35 c) 
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or 


q = T" “ + * ]T C i( S i " So) (354) 

i=i 


Fi nal ly, vising the known values of Sf(x) at d v , v = 0 , 1 , . . . , 
n + 1, one obtains the n + 1, simultaneous equations 


2 it 


(S v " s Q ) 


--h f 


Sw( d 


v ,e)dB - | 


d vN /L Q 2 -d v 2 


+ L 0 cos 


■ 1 /-dy\ 

VL 0 y. 


L n - dvd 


1 V cJcdi - d v ) 2 cosh- 1 

2 L 1 l L 0 |d v -di| 


- - Jl q 2 - d-j 2 dqcos -1 
5 L \Lo/ 


n/^q 2 - < 3 - V 2 li v = 1>2, • • n 


+ 1 


(36) 


which deter mi ne the n + 1 constants B, Cx, C 2 , . . . , Cn- These, in 
turn, fix the shape, through equation (3*0, and the wave drag, 4 through 
equations (35) . 


WINGS CENTRALLY MOUNTED ON SLENDER CLOSED BODIES OF REVOLUTION 


In the preceding section the interference between the central portion 
of the airplane and Its nose and tail regions was neglected. In this por- 
tion we will consider the entire fuselage, assuming, first, it is a slender 
closed body and, second, it can be calculated in the presence of the wing, 
using the same postulates given in the previous section under "Basic Theory 
and Assumptions 

It has come to the author's attention, through a private communica- 
tion, that W. T. Lord has obtained a solution similar to the above arirl is 
using it to calculate the drag of bodies of revolution having their ar eas 
specified at a given number of stations. Such a method has the advantage 
of giving the lower bound to the drag of bodies whose areas have been 
measured at a discrete number of places and, further, of giving a value 
representative of all area variations in the vicinity of the calculated 
optimum. 
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Unlimi ted Indentation Length, Fixed Volume 


Let us first consider the question: 

(iv) Given the wing, body length, and total volume of the com- 
bination, what is the area distribution of the body which yields 
a minimum wave drag if the apexes of the Mach forecone and 
aftercone enclosing the wing lie within the body (see sketch (f)) 
and the specified volume is large enough for the body to be real? 

This problem can be solved in a simple manner by means of 
the following le mma discovered by R. T- Jones, using methods 
similar to those introduced in reference 10 - 

Designate the closed body which has a minimum wave drag 
for a fixed volume and length as a Sears -Haack body. Then the 
total drag of a Sears-Haack body and any other wing or (slender) 
body entirely within the fore and after Mach cones with apexes 
at the tail and nose of the Sears-Haack body, respectively, is 
given by the equation 

D - %H (j + ^ ) + <37) 


drag of Sears-Haack body alone 
drag of second body alone 
V SH volume of Sears-Haack body 
Vi volume of second body 

A proof of this lemma can be obtained by placing the Sears-Haack 
source distribution and the wing equivalent multipole distributions (or 
the second body’s equivalent multipole distribution) in equation ( 6 ). 
Since only the sources interfere, the drag can be written in the form 

n P L ° r l '° a o 

D = Dgg + — / a Q (xjjdxi / + Bi (38) 

itU 0 2 01 0 Xi - x 2 

-Do _ ‘"O 

where the interference term has been integrated by parts and -L 0 T , L 0 , 
and -l c , l 0 form bounds of the arbitrary and Sears-Haack source 



Sketch (f) 
where 


I’SH wave 
Di wave 
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distributions , 


a o 1 811(1 a Ogg^ 


respectively. 


As is well known 


^SE _ ^SH 
<1 " ^o 4 


and 


_1_ 

Uo 



^SH ^ 

«1 0 4 Jl Q 2 - x 2 


(39) 


(40) 


Placing equation (4o) in ( 38 ) and integrating, one finds 


D - D sh - 


8q7g£_2_ 

V SH 



a 0i (xi) 

dxi + Di 

Uo • 


and since one can easily show 


Vi = - — / xia 0 (xi)dxi 

u o O 

~ h o 

equation (37) follows immediately. 

Returning now to problem (iv), we see that its solution follows frcaa 
equation (37) and the solution is, in fact, simply a Sears -Haack body 
having, at the appropriate place relative to the wing -body juncture , the 
additional area variation specified by equation (4). This follows, since, 
if Di represents the combined drag of the wing and indentation, then V x , 
the combined volume of the wing and indentation, is zero. Hence, the 
minimum value of D, for a given volume, is obtained when Dgg and Di are 
independently minimized. But Dgg; is already a minimum on a volume basis 
and Di is a minimum for a given wing, notice the location of the wing 
along the body is immaterial, provided the required indentation can be 
accommodated by the fuselage. 
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Limited Indentation Length On Sears-Haack Body, 

Fixed Volume 

Consider, next, the more difficult problem 

(v) Given a wing and Sears-Eaack body of length 2Z 0 (long enough 
to contain the apexes of the fore and after Mach cones enclosing 
the wing), what modification of this fuselage within the length 
+ ^1 (and within that length only, see sketch (g)) minimizes 
the total wave drag for a given total volume? 


Basic Sears - 
Haack body 



In order to answer this ques- 
tion, it is necessary to consider 
separately two cases; namely, the 
one in which Zi* > Lo' and l x > L 0 
(i.e., the portion of the body free 
for variation contains the apexes of 
the wing's Mach cone envelope, as 
shown in sketch (f)) and the other 
in which the preceding conditions 
are not satisfied. 

First consider the combination 
for which Zi * > L 0 ’ and l x >L Q . 

The wave drag of such a combination 
can always be calculated using equa- 
tion (37) wherein Dgg is the wave 
drag of the basic Sears-Haack body 
Cancellation fixed by the stationary nose and tail 
portions , Di is the combined wave 
drag of the wing and the (as yet 
u nkn own) body indentation, and Vi 
is the net difference in volume 
between V SH> the volume of the Sears- 
Haack body, and the final volume of 
the complete configuration. Since 
the basic Sears-Haack body is fixed 
Sketch (g) and the total volume is given, the 

entire term Dgg[l + (2Vi/Vgg) ] is fixed and the solution to the problem 
is obviously that for which the wing equivalent sources and the source 
simulating the body indentation combine to form a Sears-Haack distribu- 
tion in the interval -l± < x < l x . 

Using equations (2) and ( 3 ) to relate the wing and body source varia- 
tions to their respective areas, we find the fuselage cross-sectional area 
can be written for 


■Z 0 < x < -Zi 


Sf(x) 


^SH 

3*z 0 4 


Oo 2 -* 2 ) 


3/2 


(it-la) 



MCA TU 3667 


25 


for -Z x ' < x < Z ; 


8 V g 

3jtZ 0 


2Jt 


SfW = ^SIL ( lo 2 . . A. J s^WeJae 


128Vi r/ , . w, x-,3/2 

[(Zl' + x)(Zi - X) J 


3it(Z + Zi‘)‘ 


(to) 


and for Z x < x < Z D 


M*) W, 2 -* a > 3/: 


3*Z 


(to) 


The total wave drag of the wing and the fuselage, as given by equa- 
tions (4l) is then 

D = 8^ / 2VA 128 V X 2 

1 *Z 0 4 V W n(Zi+Zi‘) 4 q 

where D e is defined by equations (5)j ( 6 ), and ( 7 ). 

Since, as we have been assuming, pA is small, D e is negligible, 
and a comparison between equations ( 37 ) and (b2) shows that the drag of 
the combination formed by mounting two wing panels on a Sears -Haack body 
can be reduced without a change in the total volume and with a modifica- 
tion limited to the interval -Z x ' < x < Z x by the difference between 
the drag of the two panels flying alone and the drag of a Sears -Haack 
• body having a length equal to Z x * + Z x , and a volume equal to the volume 
of the two panels. So long as the points x = -Z x ’ and x = Z x do not 
lie off the basic body, and so long as the required indentation can be 
accommodated, this result is independent of the wing ' s fore-and-aft 
position. 

If the body modification is limited so that either Z x '< L 0 ' or 
Z x < Lo (see sketch (g)) or both, the above results do not apply, since, 
in such cases, the second body - in the sense defined above - cannot be 
varied for x between -Z x ’ and -L 0 ' or Lo and Z x or both, and its 
drag cannot, therefore, be reduced to that of an equivalent Sears-Haack 
body. The best modification in this case can be calculated from the 
results presented in the material immediately following. 
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Limited Indentation On Arbitrary Body - 
Fixed Volume 


Consider the question 

(vi) Given a wing, a body length, and the area distribution of 
the fore-and-aft portions of a body, what is the variation of 
area along the intermediate portion of the body which yields a 
mi ni mum wave drag for a fixed total volume? . 

Again, as in equation (17), let cr(x) represent the sum of the sources 
representing the basic body and the wing equivalent source distribution. 


2 jt 

r(x) = Sf (x) + ^ J Sw(x,e)d e 


m 


It is now convenient, however, to let cr(x) be a fixed function in the 
entire interval -Z 0 ’ < x < Z Q , see sketch (h) , and let the body modifica- 
tions, which are to be optimized in 
the interval -lx < x < lx be repre- 
sented by ASf(x) which has the end 
conditions 



AS f '(-lx) = ASf '(Zi) = 0 

ASpC-Zj =£S f (lx) = 0 J 




The change in volume caused by the 
body modification, AV, is given by 


AV 


■'l 

= - J' xAS-p ' (x)dx (^ 5 ) 
-Z-, 
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The usual variational procedure leads directly to the integral 
equation 



ASf" (x 2 )dx 2 
xi -x 2 



a" (x 2 )dx 2 

+ Ao + AiX,; -l x < x < 1-, 

xi - x 2 


(46) 


where and 7\ x , are fixed by the conditions given in equations (44) 
and (45). Equation (46) is similar in form to equation (13) and its 
inversion can be obtained by use of methods similar to those for invert- 
ing the latter equation. Thus, the solution to equation (46) becomes 
for -l 1 < x < l x 


AS f '(x) 


-ff’(x) 


+ ~ x 2 

it 



-tr'(xi)dxi 

“ - 1 + 

(xi - x) n/ Xi 2 - l ! 2 


I 


o 1 (xx)dxx + 4x 

l x (xi - x)*V Xi 2 - Zi 2 U 4 


2I 2 - l x 2 I 0 - 2(v + AV) 


. 2Ii 
+ lx* 


(W 


where 


In 



xi n CT'(xi)dxi 
■Jxj 2 - l x 2 


xx 11 ^ * (xx)dxx 



•J xx 2 - l x z 


m 


V is the total volume of the wing and unmodified fuselage, that is 



xcr 1 (x)dx 




and 
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Equation (V 7 ) integrates to give 

^Sf(x) = - a(x) + (l x x + li 2 Io) si Zi 2 - x 2 - 

Till 


4 

3*1 1 4 


2I 2 - l i 2 I 0 + 2(V +AV) 


( lx 2 


3/2 

X 2 ) + l(x) 

(50) 


where 


H(x) 


1 



XiX - Ii a 

7 (ll Z - x 2 )(x 1 2 - lx 2 ) 



dxi - 


1 

* 


f ff'(xi) 

r + .™-i * * 1 

J 

^ 7(^i 2 -x 2 )(x 1 2 -l 1 2 ) ^ J 


dxi 


(51) 


If Dq- is the drag of the original wing -body combination and D^o 
is the drag of a body of revolution having the same normal area distribu- 
tion as the modification, then 




+ 


8av r 

rtli 4 


2(V + AV) 


+ ^i 2 Iq ~ 



( 52 a) 


On the other hand can be written 


fire 

q. 



h 2 ! c 


- 2(V + AV) 


r--} 


( 52 b) 


where if 


G(xi,x 2 ) 


1 /x 2 2 - Ii 2 

Xl- x 2 7 x x 2 - H 2 


(53) 
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B*/q 


can be expressed as 


= 1_ 
q. 2^ 




<3^20'' (x 1 )cr tt (x 1 )G(xi,x 2 ) 


2 



dXi 



dx2tf , (xi)cf ,, (x 2 )G(xi,x 2 ) + 



dx 2 d ' (xi)a n ( x 2 )g(xi,x 2 ) 




Limited Indentation On Arbitrary Body - 
Fixed Diameter 


As a final example in tbis section, consider the question 

(vii) Given a wing, a body length, and the area distribution of 
the fore-and-aft portions of a body, what is the intermediate 
variation of fuselage area that has a given area at some inter- 
mediate station xi = di and yields a minimum wave drag for 
the combination? 

Using the same definition for cr(x) as is given in equation (43), and 
again designating the area modification as ASf(x), one can apply the 
same methods used to develop equations (2l) and (46) and write the inte- 
gral equation for. ASp(x) in the form 

h 

P AS f ,, (x 2 )dx 2 
-l! Xl - X2 




ff" (x 2 )dx 2 


+ 


It Xi - x 2 
L o 


?Vo, -lx < x < di 
Ai, di < x < lx 


(56) 


where 7^ and 7\ 1 are constants whose values depend upon the restraints. 
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The solution to equation (56) can he written 

-lx 

ASf"(x) = -a" (x) ^ [" f a - ---- Jx 2 2 - lx 2 dx s - 

-x 2 L , x -x 2 


r ° fo 2 ^ - Zi 2 dx 2 - A - Bx - Cicosh -1 — 

j x - X 2 J Zx 


Zi - xdi 


lx - dll 


(57) 


And the three constants A , B, and Cx are fixed by the conditions: 
(l) continuous slope 


J AS f "(x)dx = 0 


(2) the body area at x = Zx is unchanged 

lx 


J ' ASf *(x)dx = 0 


-lx 

and (3) the body area at dx is given 

cLi 


J ASf'(x)dx = -ASg (dx) 
-lx 


The final solution is 

£Sf(x) = -d(x) + H(x) - ~ (x - dx) 2 cosh -1 


Zx 2 - xdi 
lx lx- dxl 


(56a) 


(58b) 


(58c) 


l x 2 -dx 2 t 31 n/Zx 2 - 


Z x 2 - dxx 


dx" 


ASf(dx) + 
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where 


Ci = 


21 - 


(l i 2 -di 2 ) 2 


Jl 2 _ ^ & 

AS f (di) + cr(di) - H(di) - ■ - - -- -- (Iidi + I 0 li 2 ) 


Ttl l 2 


(60) 


and In and H(x) are defined in equations (kQ) and (51) > respectively. 
The drag can he expressed either as 


D Dft Da q 

- = JtCj/SSf (di) (6la) 

where D ff is again the drag of the original unmodified combination and 
D /V g is the drag of the modification alone , or as 



B 2 (Z x 2 -d x 2 ) 

1 + 



(6lb) 


where D* and Ci are defined by equations ( 5*0 and (60) , respectively, 
and B is given by 


B 




(62) 


Ames Aeronautical. Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 4 , 1955 
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